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Abstract 

We study the wave functions of light and heavy mesons in both hard-wall (HW) and soft-wall 
(SW) holographic models which use AdS/CFT correspondence. In the case of massless constituents, 
the asymptotic behaviors of the electromagnetic form factor, the distribution amplitudes, and the 
decay constants for the two models are the same, if the relation between the dilaton scale parameter 
and the size of meson is an inverse proportion. On the other hand, by introducing a quark mass 
dependence in the wave function, the differences of the distribution amplitudes between the two 
models are obvious. In addition, for the SW model, the dependences of the decay constants of 
meson on the dilaton scale parameter k differ; especially /q 9 ~ K 3 /mg is consistent with the 

1/2 

prediction of the heavy quark effective theory if k ~ nig . Thus the parameters of the two models 
are fit by the decay constants of the distinct mesons; the distribution amplitudes and the ^-moments 
are calculated and compared. 
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I. INTRODUCTION 



The internal structure of the hadron has continually been a very important and interesting 
subject in the research field of high-energy physics. From the perspective of the experiment, 
in terms of the measurements of decay rate or cross section, one can realize the internal 
structure; on the other hand, from the viewpoint of theory, the internal structure is displayed 
from the distribution amplitude or the form factor. No matter whether decay or collision 
happens in the hadron, the strong interaction plays an important role. Until now, the most 
successful theory for describing the strong interaction is quantum chromodynamics (QCD). 
However, it is still difficult to carry out accurate calculations based on QCD; the main key 
lies in the strong coupling constant, especially in the low-energy region, close to unity. Such 
large value makes the perturbative calculation an arduous process. 

There were several nonperturbative approaches inspired by QCD which could obtain 
some properties of the hadron and recently, based on the correspondence of string theory 
in anti-de Sitter (AdS) space and conformal field theory (CFT) in physical space-time (l- 
4j, a semiclassical approximation to QCD, light-front holography (L 
developed for describing the phenomenology of hadronic properties 
mapping of the string modes in AdS fifth dimension z to the light-front wave function in 
the impact variable (, which measures the separation of the constituents inside a hadron. 
The mapping is proceeded by matching certain matrix elements presented within the string 
theory in AdS space and the light-front theory in Minkowski space. This approach, known 
as bottom-up, allows to built models that have been successful in various QCD applications, 



7 H), was successfully 
Jl2]. LFH provides 



such as hadronic scattering processes 
and chiral symmetry breaking 



13Nl5| , hadronic spectrum 



2CH22| . and quark potentials 23N25] 
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4l9||. hadronic couplings 



It is well known that the main feature of QCD as the fundamental theory of the strong 
interaction lies in its nonperturbative behavior. The picture of nonperturbative QCD leads 
to quark confinement and chiral symmetry breaking, which are believed to be the two es- 
sential mechanisms for forming hadron formation from QCD. Undoubtedly, the complicated 
and nontrivial vacuum is not only the most important object to be understood for the non- 
perturbative aspect of QCD, but also a starting point for the construction of hadronic wave 
functions. Furthermore, in the light- front (LF) coordinates, the nontrivial vacuum is con- 
nected to the particles with zero-longitudinal momentum (the so-called zero-mode particles) 
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and the infrared divergences from the small longitudinal momentum 26] . In fact, the LF 

n 

infrared divergences are to be expected as the sources of chiral symmetry breaking [271 ] . 
However, the understanding and application of the true QCD vacuum are still very limited 
in the LF framework. Thus, here we study the relevant subjects in the LF framework with 
a trivial vacuum. An alternate approach which can reveal the phenomena of chiral sym- 
metry breaking was proposed by the authors of Ref. 2JJ, who chose the field content of 
the five-dimensional theory to holographically reproduce the dynamics of chiral symmetry 
breaking in QCD. While this holographic model depends on only three free parameters, it 
agrees surprisingly well with the seven observables of the light meson. 

For compatibility with QCD, the conformal invariance must be broken and the confine- 
ment in the infrared region must be introduced. There are two types of the modified AdS 
geometry which can achieve this results, as in the literature. One is the "hard-wall" (HW) 



model 



10. 



13 



20 



21 



28 



where partons are free inside the hadron and a hard cutoff is 



applied as the boundary. Unlike the standard bag model [31] , these boundary conditions are 



imposed on t 
(SW) model . 



re impact variable C, not on the bag of the radius. The other is the "soft-wall" 
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22 



24 
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44j which has no sharp boundaries and employs a background 



dilaton field in the AdS space as a smooth cutoff. A problem with the HW model is that 
the dependence of hadron masses on the higher orbital angular momenta is linear, which is 
different from the quadratic behavior or the so-called "Regge trajectory". Conversely, the 
SW model was initiated for solving the problem of the hadronic mass spectrum. Each of 
these two models has certain advantages which are explained the above references. 

In this work, we study the wave functions of light and heavy mesons and make the compar- 
isons of the asymptotic behaviors of the electromagnetic (EM) form factor, the distribution 
amplitudes, and the ^-moments in both the HW and SW holographic models. The case of 
the massless constituents is first considered, and a generalization regarding the behavior of 
the massive quarks 45[ follows. There are two points worth mentioning: First, for the HW 
model, we use the parameter I instead of 1/Aqcd to represent the size of the meson in both 
cases of the massless and massive constituents. Second, for the SW model, although both 



dilaton field profiles exp(±n 2 z 2 ) have been used in the literature 16|, |46J] to reproduce the 
behavior of Regge trajectories for higher spin states, and the authors of Ref. 47| related 
these solutions by a canonical transformation, the main motivation for using the dilaton field 
profiles exp(+K?z 2 ) is its seemingly better confinement properties 24|. However, as shown 
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in Refs. lg, |48[, on one hand, the dilaton exponential exp(+n 2 z 2 ) leads to the existence of 
a spurious massless scalar mode in the model; on the other hand, the desired confinement 
properties can be realized in the model with the dilaton exponential exp(— k, 2 z 2 ). Thus, here 
we apply the latter one to both the hadronic field and the bulk-to-boundary propagator. 

The remained of this paper is organized as follows. Section II reviews the extraction of 
wave functions by the holographic mapping of the light-front formulism to AdS string mode. 
In Sec. Ill, we compare both models in regard to the massless and massive constituents. In 
Sec. IV, we fit the dilaton scale parameter and the radius of bag with the decay constants of 
the distinct mesons and estimate the distribution amplitudes. The first four ^-moments in 
both models are compared with those of the other theoretical calculations. Conclusions are 
presented in Sec. V. Other useful derivations for the discussion are given in the appendices. 

II. FORMULISM 

A. AdS wave equations 

The action for a spin- J field $j in AdSs space-time in presence of a dilaton background 



field <f(z) is given by [49 1 



= j d\dz^-ge-^\d N ^jd N ^j - ti(z)<S> 2 j) } (2.1) 

where n 2 (z) = fi 2 + gj(p(z) is a "dressed" five-dimensional mass because of the interaction 
of (p(z) with $j. The mass fij is related to the conformal dimension A by (fijR) 2 = 
(A — J)(A — 4 + J) and the coupling gj will be fixed later. In addition, g = |det<7jvw|> 9mn 
is the metric tensor with M, A = 0, 1, 2, 3, 4, and the AdS metric is defined as 



d2 

ds 2 = g MN dx M dx N = —(r j ^dx fi dx u - dz 2 ), (2.2) 

z 

with R the horizon radius and rj^ =diag(l, — 1, — 1, — 1). $j is a rank tensor field 
§{x m )m 1 -m j which is totally symmetrical in all of its indices. Factoring out the plane 
wave along the Poincare coordinates we get: 

* P {x,z)j = e- ip -**{z)j, (2.3) 
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with four- momentum P M and invariant hadronic mass P M P M = M 2 , and taking the variation 
of (12.11) . the AdS wave equation for the spin- J field is: 

2" 



Z 



3-2J 



-ip(z) 



-<p( z ) 



3-2J ^ z 



d z ) + 



Hj{z)R 



$(z)j = M 2 $(z)j. 



(2.4) 



We may assume that the hadron field $p(x, z)j minimally couples to a massless vector field 
A M (x,z) with the action 5(3, 51]: 



,s'., = / d 4 xdz^e-^ z \-^F MN F MN + ig 5 A M ^ pi (x,z)jd M ^ P (x,z)j} 



(2.5) 



where F MN = d M A — d N A M and g§ is a five-dimensional effective coupling constant. In 
the AdS space, the propagation of an EM probe can polarize along Minkowski coordinates 



as: 



AJx,z) = e„e- i * x V(Q 2 ,z), A z = 0, 



(2.6) 



where Q 2 = —q 2 > and V(Q 2 ,z) is the bulk-to-boundary propagator with V(0, z) = 
V(Q 2 ,0) = 1, since we are normalizing the bulk solutions to the total charge operator, and 
a boundary limit of the external current A^(x, 0) = e fl e~ iq ' x . Taking the variation of the first 
term in Eq. ( 12. 51) . the AdS wave equation for the external current is: 



d z ) + z 2 Q 2 



V(Q 2 ,z)=0. 



(2.7) 



B. Light-front framework 



The Lorentz-invariant Hamiltonian equation for a relativistic bound-state hadron is given 



by 



P IM P^(P))=M 2 \ij(P)). j 



(2.8) 



where P M is the four-momentum in four-dimensional Minkowski space. In the LF QCD, the 
hadron four-momentum P^ can be expressed as: P = (P + ,P~,P ± ) with P ± = P° ± P 3 
and Pj_ = (P^P 2 ). The hadronic state, which is an eigenstate of P + , Pj_ and the total 
longitudinal spin J z , is normalized as 



(^(P + , P ± , J z )mP' + , Pl, J' z )) = 2(2nfP + 5. h ,j,5(P + - P> + )5 2 (P ± - P' ± ). (2.9) 
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In terms of the LF relative momentum variables Xi and k_u, the on-shell partonic momentum 
Pi can be expressed as 

pf = Xi P+, p l± = x t P ± + k ±l , (2.10) 
and the hadronic state is expanded as 

W p + , Pi ,j,)> = En/ l^d^id-E^s'iEk 



dxjCpk i 

x ^(ar^kjLi, K)\n;pi, A;), (2.11) 



where the delta functions are required from the momentum conservation, |n) is the multi- 
particle Fock eigenstates with n the number of patrons in a given Fock state, and Aj is the 
projection of the constituent's spin along z direction. The LF wave function ip n (xi, kj_j, Aj) 
complies with the angular momentum sum rules [52| J z = Y^=i + Y^i=i Li, z and is nor- 
malized as 

n r A<r a 2 \c n n 

EII / 2f2^3 ^ (167r3) ^ (1 - E^)^ 2 (E k ^)l^(x„k ±t , A,)| 2 = 1, (2.12) 



n,Ai 8=1 v j=l i=l 

which is derived from (12. 9p . 

In the LF framework, a lot of physical quantities relate to the LF wave function. Here we 
show some of them which are of concern. First, the decay constant of a pseudoscalar meson 
f M is given by: 

(0\J+\M(P + ,P ± ))=iP + f M , (2.13) 

where = gi7 + (l — 7 5 )q , 2 is the flavor changing weak current which is evaluated at fixed 
light-cone time x + = 0. Only the valence quarks contribute to the decay; thus, we expand 
the initial state in Eq. ( 12. 13f) into the Fock component and find 

f M = 2^2N C J ^^ftCn, k x ). (2.14) 

Equation (12.14p enables the straight evaluation of the decay constant in terms of the LF 
wave function ipq q (x, kx). 

Next, the distribution amplitude <f>(x, Q), which is the amplitude for finding constituents 
with longitudinal momentum fraction x in the meson and collinear up to the scale Q, is 
defined as 



53j: 



M)= I |^Mx,k ± ). (2.15) 



The meson state \M), which is expanded into Fock states, will approximate \qq) when Q 2 is 
large. Therefore, we can also calculate <fr{x) = <p(x, Q — > oo) by the wave function ipq q (x, kj_). 
In addition, an approach to parameterize distribution amplitude is to calculate the so-called 
^-moments: 

n = ri JL Z > (2 - 16) 

J_i <w0 

where £ = 1 — 2a;. 

Finally, the EM form factor of a meson is defined as: 

(M(P')\J + (0)\M(P)) = (P + P') + F(Q 2 ), (2.17) 

where J + (y) = J2 g e q ( l(y)l +( l(y)^ P' = P + q, and F(0) = 1. Using the expanded state Eq. 
( 12. lip and the normalization condition Eq. f)2.9p . the EM form factor of a meson can be 



expressed as [54-1561] : 



F(Q ) = / dxp(x,q ± 
dxd 2 k 



2(2tt) 



^e i ?/'i (? (x,k_ L )V'g (? (x,k_ L ), (2.18) 



where p(x, q ± ) is an effective single-particle density and kj_ = kj_ + (1 — a?j)q_i_. Equation 
(I2.18P reveals that the current matrix element in Eq.f l2.17P can be represented as overlaps 
of the LF wave functions. 



C. Light-front holographic mapping 

For the holographic mapping of the light-front wave function to AdS string mode, it is 
convenient to define the transverse center of momentum of a hadron as: 

R± = j dy-d 2 y ± T ++ y ± , (2.19) 

where is the energy momentum tensor. Then the partonic transverse position r±i and 
the internal coordinates bj_j, which conjugate to the relative momentum variable kj_j, have 
the following relations: 

n n 

XiV ±i = XiR± + b_Li, 2j b -Li = °> Rl = /"l XiT±i. (2.20) 

i=l i=l 
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The light-front wave function tp n {xi, k_u) can be expressed by the internal coordinates b_u 
as: 

n— 1 „ n— 1 

^n(^,k ±l ) = (4tt)^ 1 )/ 2 n / d 2 b ±4 exp(^b ±J • k ±l )^ n ( Xl ,h ±l ), (2.21) 

where ip n {xi, b±i) is the light-front wave function in the coordinate space and is normalized 
as: 

n „ 

J2U / dx i( i 2 b ±i \ip n (x h b ± i)\ 2 = 1- (2.22) 

n i=l 

Then, we can substitute Eq. f)2.2ip into Eq. (I2.18P and integrate over kj_ phase space. The 
EM form factor of a meson F(Q 2 ) can be obtained as: 

F(Q 2 ) = J dxd 2 b ± exp(iq ± • (1 - x)b ± )\^ qq (x, h ± )\ 2 . (2.23) 

We can also express F(Q 2 ) in terms of an effective single-particle transverse distribution 
p(x,c ± ) I 



F(Q 2 ) = J dxd 2 c ± exp(-ic ± ■ q ± )p(x,c±), (2.24) 

where Cj_ = (1 — x)bj_ is the x- weighted transverse position coordinate of the spectator 
quark, and 

p(x,c ± ) = J ^pexp(ic ± ■ q ± )p(x, qj 

= J d 2 b ± 5 2 ((l-x)b ± -c ± )\^ q (x,b ± )\ 2 - (2.25) 
If we integrate Eq. (12.241) over angles, the form factor can be obtained as [8|]: 



F(Q 2 ) = 2n j dx 1 —^ j CdCjAcq\J l —^\p{xA), (2.26) 



where ( = x{\ — ~x)\b±\ and J n is the Bessel function of the first kind. 

On the other hand, if we consider a pseudoscalar meson coupled to an external EM field 
in AdS space, the second term of Eq. (12. 5 1) can be related to a hadron matrix element as 

OS 

J d^xdz^e'^A M ^ pi {x, z)j,d M § P (x, z)j 

~ (2vr) 4 5 4 (P' - P - q)eJM(P', J')\J»\M(P, J)). (2.27) 



Substituting Eqs. (I2.3P and ( 12. 6 j) to the left hand side of Eq. ( I2.27P and extracting a delta 
function from momentum conservation at the vertex, the form factor can be obtained as: 

F(Q 2 ) = R 3 J ^L e ^$(z)V(Q 2 , z)®{z), (2.28) 

and be treated as the overlap in the fifth dimension coordinate z of the normalizable 
modes which correspond to the incoming and outgoing mesons, $p and <3>p/, with the non- 
normalizable mode, V(Q 2 , z), which deal to the external source. Comparing Eq. ( 12. 28ft with 
Eq. (12.26P and identifying the holographic variable z with the transverse LF variable (, the 
holographic mapping of the light-front wave function to AdS string mode can be accom- 
plished. Next, two kinds of dilaton background fields were considered and their respective 
wave functions extracted. 



1. Hard Model 



The simplest dilaton background field is none; that is, <p{z) = 0. Thus, the AdS wave 
equation, Eq. (12.71) . for the external current is reduced as: 

[z 2 d 2 z -zd z -z 2 Q 2 }V H (Q 2 ,z) = 0, 



and its solution with the boundary condition is: 

V R (Q 2 ,z) = zQK 1 (zQ), 



(2.29) 



(2.30) 



where K a (x) is the modified Bessel functions of the second kind. In addition, an integral is 
done I8f as: 



dxJ CQ 



1 — x 



X 



(2.31) 



Substituting Eqs. (I2.30p and (I2.3ip to Eqs. (12.281) and (12.261) . respectively, and comparing 
the latter two equations, the relation between the LF wave function and the AdS string 
mode is obtained as: 



1 7 t ni^ ^ n J $ h(C)I 2 



(2.32) 



2tt"^ -' C 4 

where z = ( and Eq. (I2.25P are applied. On the other hand, for the meson field, Eq. (12. 4p 
is reduced as: 



2 3 - 2 J / \ijR 
-o z H d z 4 



$ H (*) = M 2 $u(z). 



(2.33) 



By substituting $h(C) = (C/R) 3 ^ 2+J( &h(C); an effective Schrodinger equation is obtained 
as: 

1" 



d 2 4L 2 
+ 



$h(C) = -M^h(C), 



(2.34) 



d( 2 4C 2 

where the conformal dimension A = L + 2 and (fijR) 2 = L 2 — (2 — J) 2 are used. However, for 
the known confinement inside a hadron, a baglike model 3l| is considered where partons are 
free inside the meson and are forbidden outside the meson. An additional hard-well potential 
is needed: U(() = if ( < I and U(Q — oo if £ > I, where I is the size of hadron. Then, the 
solution to Eq. (12. 34p is $h(C) = C( 1/2 J L ((M) with a boundary condition <f> H (C = I) = 
and a normalization condition J d(\&(()\ 2 = 1. As for the meson field, the solution to Eq. 
(12331) is: 

z 2 ~ J V2 



-Jl(zM). 



(2.35) 



R 3 / 2 - J U L+l {Ml) 

From Eq. (I2.32p . the LF wave function in the HW model in the limit of massless constituents 
is: 

1 r-. r . . r~, r„ . , . /. , I 2 



1pL,n(x,bj 



\fTtlJ L+ x{P Lyn 



V*(l - x)J L (y/x(l - x)\b ± \M L , n )6 b 2 ± < 



X 1 



X 



(2.36) 



where (3L, n is the root of the Bessel function, M.L,n = &L,n/h an d 6 is the step function. 



2. Soft Model 



Another background field which introduced an infrared soft cutoff is a dilaton ^p(z) 
k 2 z 2 . This time Eq. (12. 7\i is reduced as: 



z 2 d 2 z - (1 + 2n 2 z 2 )zd z - z 2 Q 2 }V s (Q 2 , z) = 



and the solution with the boundary condition is: 



^ 2 ^) = r(i + ^)t/(g,o,.V 



(2.37) 



(2.38) 



where T(x) is the Gamma function and U(a,b,c) is the confluent hvpergeometric function. 
In the large Q 2 limit, that is, Q 2 ^> 4k 2 , the solution is reduced as 8|: 

V S (Q 2 , z) zQK^zQ) = V u (Q 2 ,z). (2.39) 
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Thus, in the large Q limit, the relation between the LF wave function and the AdS string 
mode is obtained as: 



R 3 

For this background field, Eq. f!2.4[) is reduced as: 



-kvI^sCOI 1 



C 4 



(2.40) 



-d 2 + 



3-2J + 2k 2 z 2 ffij(z)R 
<9 2 + 



(2.41) 



and an effective Schrodinger equation can be obtained as: 

rP A T 2 — 1 

+ ^—± + K 4 z 2 + K 2 {gjR 2 -2(J-l)} 



$s(C) = A4 2 $s(C), (2-42) 



d( 2 AC 2 

where $ S (C) = e- K ^ 2 / 2 (C/i?)" 3/2+J $s(0- The normalized solution to Eq. (I2^2]l is: 



MO = « 



2n! 



(n + L)! 



e^^L^C 2 ) 



(2.43) 



with 



A .2 „ 2 / ^ L- J + 2 + gj R 2 /2 
A4 = 4k n H 



(2.44) 



and the meson field is: 



2 



2-J+L 



i? 3 /2- 



-K 



i+1 



2n! 



(n + L)! 



i*V). 



(2.45) 



In order to obtain a massless pion and a linear Regge trajectories: Ai 2 ~ J at large J, the 
value gj-R 2 = 4(J— 1) is fixed. Thus, from Eq. (12.401) . the LF wave function in the soft-wall 
model in the limits of massless constituents is: 



V> L (x,bj 
with 



K 



L+l 



111 



[x(l - x)] {L+1)/2 \b ± \ L e~ K2x(1 ~ x)h ^ /2 L^K 2 x(l - x)b 2 L ),(2.46) 



M 2 = Ak 2 in + 



L + J 



(2.47) 
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III. TWO HOLOGRAPHIC MODELS IN COMPARISON 



massless constituents 



In this section, we make the comparisons between the HW and SW models for some 
hadron properties. For the EM form factor of the charged pion, we take n = J = L = and 
substitute Eqs. fl2T30|) . f l2~35]) and f l2~38]) . ( |2^5|) into Eq. f l2~28j) : 



2 / dz- 



T 2 I Z Po,l 



zQK^zQ), 



f!(Q 2 ) = 2 J to 2 (L°(.V)) 2 r {i + ^)u(^ 2 , o, «v) , (3.1) 

which correspond to the HW and SW models, respectively. For the HW model, the analytical 



result is not easily obtained. As for the SW model, the EM form factor is 

Ak 2 



8]: 



Ak 2 + Q 2 ' 

However, if we check the asymptotic behavior in the large Q limit, the results are 

Qu 



(3.2) 



F n (Q 2 ^oo)^ 



R 3 e^$(z) 2 



Q 1 



(3.3) 



where u = Z(oo) for the HW(SW) model. Thus: 

2 



Q 2 F*{Q^ 



4 



Q 2 F*{Q' 



(3.4) 



JJi(/3 ,i) 

They are derived in Appendix A and, of course, the latter is easily checked according to Eq. 
(13. 2p . In other words, if we set a parametric relation: 

K= Lh(M' (3 ' 5) 

the asymptotic behaviors of Q 2 F 7T (Q 2 ) for both models are the same. In addition, the mean 
square radius of the meson P is determined from the slope of Fp(Q 2 ) at Q 2 = 0: 

dF P (Q 2 ) 



(r P ) 



-6- 



(3.6) 



Thus, from Eq. f)3.2p . the mean square radius of the pion in the SW model is 



2\S 



(3.7) 
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For the HW model, although the value of {r 2 ) R diverges logarithmically, as mentioned in 
Ref . [8J , this problem in defining mean square radius of the pion does not appear if one uses 
Neumann boundary conditions, and (r 2 ) n ~ I 2 . Thus, the inverse relation between k and I 
is still satisfied. 

One may speculate whether or not there are other similarities between these two models. 
Thus, we consider $h,s(C)> which are the normalized solutions of the effective Schrodinger 
equation and are derived from the meson field, for L = 0. If the parametric relation ( 13. 5 p is 
used, we have: 

oo 

$r°(c) = E c ^s(0, (3.8) 

n=0 

where |c n | 2 = 1. It should be noted that c n is independent of the parameters k and Z, 
which is described in Appendix B. Qualitatively, the coefficient c tells us "how much" $s(A) 
is contained in $^(A). We easily find cq ~ 0.980. In other words, the former is dominated 
within the latter. 

It is convenient to study other hadron properties in the momentum space. The LF wave 
function of the HW model for n = L = can be obtained from Eqs. (I2.2ip and (I2.36P as: 



//is 4tt v /x(1 -x)p 0tl l ( kl , 

tpn{x J k ± ) = — — — J — = , (3.9) 



X) 



where k = |kj_|. The distribution amplitude and decay constantcan be obtained from Eqs. 
(l2~T5j) and (l2TT2j) . respectively: 

y/x(l-x) V6 1n v 

Mx) = mmm- h = sU^T) ■ (3 - 10) 



and the ^-moments can be obtained by Eq. ( 12. 161) : 

(e )H " 20F F[f + 2]- (3 ' n) 
We see that the ^-moments is zero when iV is odd. In other words, the distribution amplitude 



is symmetric for £ = or x = |. This is apparent because of ^h(^) oc yx{\ — x). 

On the other hand, also for L = 0, the LF wave function of the SW model in momentum 
space can be obtained from Eqs. ( 12.211) and (I2.46P as: 

k ±) = - e^^Ll — — , (3.12) 

Ky/X(l-X) \X(1-X)K 2 J 
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and the distribution amplitude and decay constant for n = can be obtained as: 

y/x(l-x)K V6k 
M x ) = ~ ^ , fs = -g~- (3.13) 



Equations (I3.10p and 03/13]) are consistent with the result of Ref. js] where they used 
a derivation: (f)(x,Q — > oo) — >■ ^(a;, bj_ — )■ 0)/\/47r as £ — >■ 0. The ^-moments of this 
distribution amplitude, because of the same dependence as that of the hard-wall model: 



4>s(x) oc y/x(l — x), is (£, N )s = (£, N )k- m fact) if we again set the parametric relation in Eq. 

(BSD , then: 

0s(z) = 0h(x), / s = / H . (3.14) 
As for the similarity between V'H^jkj.) and ^>g(x, kj_), we can follow the similar argument 



for $h,s(C) i n Appendix B: introducing a new variable k' ± = kj_/ y/x(l — x)k, and rewriting 
the LF wave functions of the two models as 

^ (k/±) = ^wiMPoT-kVJi(M) Jo \*ikd) ' (3,15) 



and 



^s(k±) = "i 6 "^' (3.16) 

'71 



where k' =\ k' ± | and both new functions satisfy the normalization condition: J d 2 lt ± 
^'(k'j_) | 2 = 1. We may consider a complete set: {'?/'^ IO (k / ± )} which is the solution of the 
Schrodinger equation for a two-dimension harmonic oscillator: 

< H °(kl) = t^ e ~4 L ^). (3.17) 



/7T 

Here the function ip'sO^x) = V'o (kj_), which is the "ground state" for this two-dimension 
harmonic oscillator, and the function "^(k'jj can be written as a linear combination of 

{^ HO (k' ± )}: ^H(k' ± ) = E« H °( k ±)- We can easil y check that c o - °- 980 > c 'n = (-!) n Cn, 
and X} c n — 1- The coefficients also satisfy the relation: 

M^±) = Y, C 'n^±)- (3-18) 
/I 

This means that the SW wave function is dominated within the HW one for n = L = 0. 

We can realize this situation in terms of sketching ^g(k'j_) and ^(k'jj in Fig. [IJ If the 
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FIG. 1: Profile of the LF wave functions for the massless constituents. The solid and dotted lines 
correspond to V's(^J-) an d ^hO^-l)' respectively. 

integration of for both sides of Eq. (13.181) is taken, we obtain 



y/x(l — x) sr-^ , — x)k 



2tt 



(3.19) 



where 



/ 



(3.20) 



is independent of n. Thus, using the parametric relation Eq. (13. 5j) . we have another con- 



straint of c' n : 



(3.21) 



The convergence of c' n is slower than that of c n because the sign of c' n may be positive 
or negative. 

B. massive constituents 



A simple generalization of the LF wave function for massive quarks follows from the 
assumption that the momentum space LF wave function is a function of the invariant off- 
energy shell quantity [45]: 



M 2 -e 



x 



1 — x 



o • 



(3.22) 
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In other words, one may make the following replacement: 



+ m. 



x 1 — x 



M . ' -'-'0 T '"12- m 12 — 1" i 5 (3.23) 

x(l — x) x(l — X) x 1 — X 

where m; is the current mass of quark. For the SW model, this replacement is equivalent to 
a change of the kinetic term in the effective Schrodinger equation Eq. (I2.42p 491 ] : 



d 2 



+ m\ 2 . 



(3.24) 



The alterant masses of the light and heavy mesons were obtained by combining Eqs. ( I2.42p 
and (I3.24p and solving them: 

L + r 



Mi 



+ N I dx m? 9 e 



2 

'12 



(3.25) 



2 
•' L2 



where N is the normalization constant fixed from the integral N J dxe~~^~ = 1. 

Here we generalize this replacement to the HW and SW models. Thus, the LF wave 
functions of these two models for n = can be replaced as: 



ip s (x,k ± ,mi 



y/x(l-x) (^-Mo 2 ^) 

47T _ M o 



MM l), 



Kyf x{\ — X) 

If the distribution amplitude for massive quarks is considered, we find that they are just the 
integrations of the massless LF wave function with an infrared cutoff: 



-e 2k 



(3.26) 



oo 72 



d 2 k ± 

2(2tt) 



7 ip HtS (x,k ± ,mi) 



kdk 



2(2tt 



where /z = \J x{l — x)mi2. For the SW model, we obtain 



\/x(l — x)K - m i2 

(ps{x,mi) ~ - - e 2^ 



which is consistent with the result of 
obtain: 



2tt 



(3.27) 



(3.28) 



49J. As for the HW model, we may use Eq. (I3.18P to 



J u 



kdk 



2(2tt) 2 

Thus, the distribution amplitude for the HW model is 



C(x,k ± ). 



(3.29) 



\/x(l — x)k ~ m i2 



2tt 



c' + c ; i + 



in 



12 



/ / j , ^12 



+ cU 1 + 



m l2 m l2 , m l2 



2k 4 + 6k 6 



+ cW 1 + 



r?) 



12 



m 



12 



+ 



fit 



12 



3k 6 24^ 



+ . . . 



(3.30) 
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Because of the constant term of L^{t), for all n's, being equal to 1, we substitute Eq. ( I3.2ip 
into Eq. (I3.30P and obtain: 



!> H 0,m;) ~ <f) S (x,mi) 



m l2 , / m l2 , / / m l2 m l2 , m l2 



I + c ' 1ZI 2 , / ^12 , / ^12 _ ^12 , 
i + t l ..2 +t 2< 1 ..l +t 3 ..2 0.-4 ^ 



(3.31) 



k 4 3k 6 24k 8 , 

Thus, the difference between 4>H(x,mi) and <fis(x,mi) is displayed as a function of m\ 2 /^ 2 - 
If the ratio m 2 2 /K 2 <C 1, that is, the meson is composed of the light quarks, the predictive 
distribution amplitudes from the HW and SW models are nearly the same. In contrast, if 
m i2/ K ' 2 > I; the parameter relation k = l/ZJi(/3 ,i) may be unable to be satisfied and the 
distribution amplitudes will be quite different. We will display these comparisons in the 
next section. 

In addition, for the SW model, if we consider the heavy and light quarkonium states, the 
scalings of the leptonic decay constants are: 



fQQ T72 - ' /« 5 [a 1 + a 2 \n — ) — , (3.32) 



in 



m 



K 



which is described in Appendix C. The latter is consistent with the case of the massless quark 
in Eq. (13. 13ft : / ~ k. On the other hand, if we take the heavy quark limit (rriQ — > oo), the 
scaling of the leptonic decay constants of heavy-light mesons is: 

fQq ~ ~> ( 3 - 33 ) 
m Q 

which is described in Appendix D. This result is in agreement with the heavy quark effective 
theory (HQET): 

/hqet ~ 1/rriQ, (3.34) 

1 12 

if k oc tjiq . These results reveal that the dilaton scale parameter k seems to vary with the 
quark mass; this inference is in accordance with the conclusion in Section III. D of Ref. 49], 
even though they fixed k and varied an additional parameter \ q Q with the distinct mesons. 



IV. NUMERICAL RESULTS AND DISCUSSION 



For the pion, whose constituent quarks approximate masslessness, the numerical results of 
the EM form factor for the HW and SW models are obtained by Eq. (13.11) and the parameters 
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10 



<2 z (GeV z ) 



FIG. 2: Q 2 F n (Q 2 ) as a function of Q 2 . The solid and dotted 



models, respectively. Data are taken from 



661 ] (triangles), 



ines correspond to the SW and HW 



671 ] (circles), and [68] (boxes) for large 



Q transfers. The long dashed line is the limiting behavior 4k 2 . 



l n = 0.892 fm, n n = 0.426 GeV, which fit the decay constant of pion / H = fs = f w — 130.41 



MeV 57] . and shown in Fig. [2j We find that the difference between the two lines is large 
for the region lGeV 2 < Q 2 < 4GeV 2 , while for the others, it is relatively small. The long 
dashed line in Fig. [2] is the limiting behavior 4k, 2 for both models. In addition, the line 
for the SW model is a little different from that in Fig. 2 of [8J because for the latter, the 
parameter k = 0.375 GeV is obtained by fitting the data for the form factor. However, for 
the parameter = 0.426 GeV, the mean square radius of pion (r 2 ) s = 0.32 fm 2 , which is 
small compared with the PDG value (r 2 ) = 0.45 fm 2 [571 ] . 

For the other light and heavy mesons, we apply the values of current quark mass as 
m u = 2.5 MeV, m d = 5.0 MeV, m s = 100 MeV, m c = 1.29 GeV, and m b = 4.19 GeV 0. 
As for the dilaton scale parameter k, we try to let it change with the decay constant of the 
different mesons. Here we use the following values: fx- = 156.1 MeV, fr>+ = 206.7 MeV, 



f D + = 257.5 MeV, f Bd = 193 MeV, f Ba = 253 MeV |57j, f Vc = 335 MeV [5j, f Bc = 489 
MeV, and f Vb = 801 MeV | l59fl . From the above input values, the parameters k and I can 
be fixed as shown in Tabled! We find that, except for the ir and K mesons, the parameter 
relation k = l//Ji(/3o,i) is no longer satisfied. For the light meson, the last subsection shows 
that (r 2 ) s oc k~ 2 and (r 2 ) H oc I 2 . We compare the experimental data: (r 2 ) / ' (r 2 K+ ) ~ 1.44 
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TABLE I: Parameters k and I for the various mesons 





TT 


K 


D 


D s 


B 


5, 




Vb 


£ c 


k (GeV) 
I (fm) 


0.426 
0.892 


0.503 
0.745 


0.909 
0.337 


1.03 
0.306 


1.49 
0.143 


1.67 
0.139 


1.18 
0.222 


2.99 
0.0751 


2.08 
0.112 




0.0 0.2 0.4 0.6 0.8 1.0 



FIG. 3: Distribution amplitudes of the light mesons. The solid line corresponds to the tt meson 
for the SW and HW models. The long dashed and dotted lines correspond to the K meson for the 
SW and HW models, respectively. The latter two are nearly the same. 



571 1 with the ratio: 

~ 1.40, ^ ~ 1.43. (4.1) 
In addition, for the heavy-light meson, we compare the square root of heavy quark ratio 



^Jm c jm^ ~ 0.555 57 1 with the ratios 



— ~ 0.611, — ~ 0.618. (4.2) 

Thus, we can conclude that these parameters are consistent with the data and satisfy the 
result of HQET, Eq. f !3.34|) . Next, we use these parameters to evaluate the distribution 
amplitude of the mesons. The normalized distribution amplitudes of the light mesons (tt, K) 
for both models are plotted in Fig. [3J As shown in Eq. f)3.3ip for the light meson, the 
normalized distribution amplitudes for the SW and HW models are nearly the same. As 
for the other mesons, the normalized distribution amplitudes of the heavy-light mesons D, 
D s , B, and B s for both models are plotted in Figs. SI , O El and El respectively. The 
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0.0 0.2 0.4 0.6 0.8 1.0 

X 

FIG. 4: Distribution amplitudes of D meson. The solid and dotted lines correspond to the SW 
and HW models, respectively. 

2.5 | 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 




-0.5 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

0.0 0.2 0.4 0.6 0.8 1.0 

X 

FIG. 5: Distribution amplitudes of D s meson. The solid and dotted lines correspond to the SW 
and HW models, respectively. 

normalized distribution amplitudes of the heavy quarkonium states r] c , rjb, and B c for both 
models are plotted in Figs. E\ , El and [HJJ respectively. These figures show that the 
normalized distribution amplitudes of the heavy-light and heavy quarkonium mesons for the 
HW and SW models are discriminating. Finally, we evaluate the first four ^-moments of 
these distribution functions, (^ 7V ) S,H , as defined in Eq. (I2.16p . and compare the results with 
the other theoretical calculations in Table [Til 
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FIG. 6: Distribution amplitudes of B meson. The solid and dotted lines correspond to the SW 
and HW models, respectively. 
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0.0 0.2 0.4 0.6 0.8 1.0 

X 

FIG. 7: Distribution amplitudes of B s meson. The solid and dotted lines correspond to the SW 
and HW models, respectively. 

V. CONCLUSIONS 

We have compared two types of wave functions for pseudoscalar mesons in the light- 
front framework, obtained by the AdS/CFT correspondence within the hard- wall and soft- 
wall holographic models. In the case of massless constituents, we find that the asymptotic 
behaviors of Q 2 F W (Q' 2 ) : the distribution amplitudes, and the decay constants for both models 
are the same if a parametric relation, k = 1// Ji(/3 ,i), is set. Furthermore, in terms of the 
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FIG. 8: Distribution amplitudes of r/ c meson. The solid and dotted lines correspond to the SW 
and HW models, respectively. 
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0.0 0.2 0.4 0.6 0.8 1.0 

X 

FIG. 9: Distribution amplitudes of % meson. The solid and dotted lines correspond to the SW 
and HW models, respectively. 

normalized wave functions of the SW model as a complete set, the ground state of the SW 
wave function dominates within that of the HW one. On the other hand, by introducing 
a quark mass dependence, the differences of the distribution amplitudes between the two 
models are obvious, and the above parametric relation is no longer satisfied if the decay 
constants of the various mesons are regarded as inputs. In addition, for the SW model, 
the dependences of the decay constants of meson on the dilaton scale parameter differ: 
f qq ~ k + 0(m 2 q /K), f QQ ~ K s/2 /m^ 2 , and f Qq ~ ^/m 2 Q . The last one, if k ~ m^ 2 , is 
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X 

FIG. 10: Distribution amplitudes of B c meson. The solid and dotted lines correspond to the SW 
and HW models, respectively. 

1 /2 

consistent with HQET: fg q ~ l/ m g • Thus, we fit the values of k and I with the decay 
constants of the distinct mesons and find that the ratios of parameters are consistent with 
the prediction of HQET and with the ratios of the mean square radius for the light mesons. 
Finally, we plot the distribution amplitudes of mesons for the two models and compare the 
first four ^-moments of our estimations with those of the other theoretical calculations. 
Acknowledgements 
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Appendix A: EM form factor in large Q limit 

From Eq. (12.281) . we have the EM form factor of pion as: 

F„(Q 2 ) = R 3 ^e^$(z)V(Q 2 ,z)$(z), (Al) 
Jo z A 

where u = Z(oo), <p(z) = 0(-kV), and V(Q 2 ,z) = V H (Q 2 , z)(V s (Q 2 , z)) for the HW (SW) 
model. Recalling Eq. fl2T39|) . that is, V s (Q 2 ,z) ->■ V R (Q 2 ,z) = zQK x {zQ) in the large Q 
limit, we rewrite Eq. ( 1A1I) as: 

F «(Q ) = 75 / P-izQfK^zQWzQ), (A2) 

v Jo z 
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TABLE II: First moments of the distribution function (£ } for the various mesons 
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The integral Eq. f ]A2j) can be taken the integration by parts, and the result is: 
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FIG. 11: Meijer G function G^zQ / 2, 1/2). The dotted line equals 1. 

where ty(z) = e Lpi - z ^(z) 2 / z 4 , W(z) = dty (z) / 'd(zQ) , and G™'™(z,r) is the generalized form 
of Meijer G function. We plot the curve of this Meijer G function in Fig. [TT] and find it 
approaches 1 and in the large Q limit and Q — » 0, respectively. Thus, the first term of Eq. 
( 1A3I) vanishes because $(2) ~ for both two models, and the second term approximately 
equals a constant when the large Q limit is taken: 



large 



R 3 



*'{z)d{zQ) 



R 3 e v M$(z) 2 



10 Q 

Substituting the relevant functions and numbers, we find: 



Q 2 F?(Q 2 ^ 



4 



/Ji(/3o, 



(A4) 



(A5) 



Appendix B: Correlation between c n and k 



The solutions of the effective Schrodinger equations, Eqs. (I2.34p and (12.42 p . for L = 



arc: 



and 



*h(0 



2C j ( CA),n+l 



I 



n(0 = ^e-^ 2 /"Ll{n\% 



(Bl) 



(B2) 



where f c?CI ( ^ ) h,s(C)| 2 = 1 is satisfied. If Eq. (13. 5p is substituted to Eq. (13.81) . we have: 

'l/reJiGSo.O 

dCV^C«^o(C«^ ,iJi(^o,i)) x ^^~ kH2/2 LI{k 2 C). (B3) 
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We may introduce a variable A = k 2 ( 2 , then $h,s(C) can be rewritten as 

<(A) = Jo(VA/V+i^i(/W)), 
$' S "(A) = e-^L° n (X), 



(B4) 



which satisfy the normalization condition J <iA|$ /n (A)| 2 = 1. It is well known that the 
associated Laguerre polynomials L^{X) have an orthogonality: 



(B5) 



e- x L° n (X)L° m (X)dX = 5 m , n . 
Thus, {^"(A)} is a complete set and Eq. (1B3j) can be rewritten as: 

-l/JiWo,i) 2 

dXJ (V\p 0> iJi(M) x e- x ' 2 Ll{X). (B6) 
It is obvious that c n is only dependent of (3 and is independent of the parameters k and I. 



Appendix C: Scalings of /qq and f qq 



In considering the heavy quarkonium, the LF wave function for the SW model is: 



i/; s (x, k_L,mi) = N ( 
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where Nqq is the normalization constant and 
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If we change the variable sin 9 = 1 — 2x, the normalization constant is replaced as: 
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The distribution amplitude and the decay constant can be evaluated as: 
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where r = itlq/ 'k and erfc(a) is the complementary error function. Performing an expansion 
in powers of 1/r for Eq. flC5j) . we obtain: 



+ o 



-5/2 



(C6) 



27r 3/4 r l/2 

Taking the heavy quark limit, we have the scaling: J'qq ~ K, 3 ^ 2 /mQ 2 . For the light quarko 
nium, we perform an expansion in powers of r for Eq. ( 1C5I) and obtain: 



gg 



\/Qk k 13 
~8 4 V 2 



3 + 2 7 + 21n[2r] + T' 



+ O (r 3 ) , 



(C7) 



where 7 is the Euler's constant and T'[a] is the digamma function. 



Appendix D: Scaling of fq q 



In considering the heavy-light meson, the LF wave function for the SW model is: 
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Taking the limit mq 3> m q and ignoring the light quark mass, we obtain: 
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K K 

where F[a, z] is the incomplete gamma function. Thus, the distribution amplitude and the 
decay constant can be evaluated as: 
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We perform an expansion in powers of 1/r for Eq. (1D5|) and obtain: 



(D5) 
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Taking the heavy quark limit, we have the scaling: /q 9 ~ tc" /tuq 
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